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Abstract 

We use very deep near-infrared imaging data taken with Multi-Object InfraRed Camera and 
Spectrograph (MOIRCS) on the Subaru Telescope to investigate the number counts of Distant Red Galaxies 
(DRGs). We have observed a 4x7 arcmin 2 field in the Great Observatories Origins Deep Survey North 
(GOODS-N), and our data reach J=24.6 and K=23.2 (5<r, Vega magnitude). The surface density of DRGs 
selected by J- K > 2.3 is 2.35±0.31 arcmin -2 at K < 22 and 3.54±0.38 arcmin" 2 at K < 23, respectively. 
These values are consistent with those in the GOODS-South and FIRES. Our deep and wide data suggest 
that the number counts of DRGs turn over at K ~ 22, and the surface density of the faint DRGs with 
K > 22 is smaller than that expected from the number counts at the brighter magnitude. The result 
indicates that while there are many bright galaxies at 2 < z < 4 with the relatively old stellar population 
and/or heavy dust extinction, the number of the faint galaxies with the similar red color is relatively small. 
Different behaviors of the number counts of the DRGs and bluer galaxies with 2 < z p hot < 4 at K > 22 
suggest that the mass-dependent color distribution, where most of low-mass galaxies are blue while more 
massive galaxies tend to have redder colors, had already been established at that epoch. 
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1. Introduction 

The Distant Red Galaxies (DRGs) are selected in the 
near-infrared (NIR) wavelength by the simple criterion 
J — K > 2.3 in order to sample high-redshift (z > 2) galax- 
ies with significant fraction of evolved stars such as nor- 
mal galaxies seen in the present universe, which are of- 
ten missed by the "drop-out" selection technique for the 
Lyman break galaxies (Franx et al. 2003) . These red J — K 
colors are produced by Balmer / 4000A-break at 2 < z < 4 
(the breaks enter between J and if-bands) and/or heavy 
dust extinction at similar redshift. In fact, several DRGs 
are confirmed to be at 2 < z < 4 by spectroscopy (van 
Dokkum ct al. 2003, Rcddy et al. 2005), and their photo- 
metric redshift also lies in the range of 2 < z < 4 (Forster 
Schreiber et al. 2004). 

Several studies found that many DRGs have high star 
formation rates (van Dokkum et al. 2004, Knudsen et 
al. 2005, Reddy et al. 2005), while some of them seem 
to have little star formation activity and evolve passively 
(Labbe et al. 2005, Papovich ct al. 2006). Analyses of the 
spectral energy distributions (SEDs) suggest that DRGs 
are more massive (M* > 10 11 M Q ) and older (~l-3 Gyr 
old) than the UV-bright LBGs at similar redshifts (Forster 



Schreiber et al. 2004, Iwata et al. 2005, Labbe et al. 2005, 
Papovich et al. 2006, Kriek et al. 2006). Recently, the 
strong angular clustering of these galaxies has also been 
reported (Grazian ct al. 2006, Foucaud et al. 2006, Quadri 
et al. 2006). 

However, most of these previous studies have focused 
on the bright (massive) part of DRGs. It is important to 
sample fainter DRGs and investigate their properties in 
order to understand the formation and evolution of more 
"normal" galaxies. Galaxies with the luminosity of L* or 
sub-L* at z ~ 3 would have K-band magnitude of K ~22- 
23. Such very deep NIR data have been limited to rela- 
tively small areas so far (e.g., Maihara et al. 2001, Labbe 
et al. 2003, Minowa et al. 2005). 

In this paper, we use very deep and wide NIR data 
taken with MOIRCS on the Subaru Telescope in order 
to investigate the number counts of DRGs down to K = 
23. The wide field of view of MOIRCS (4x7 arcmin 2 ) 
and the large collecting area of the telescope allow us to 
analyze DRGs to the faint magnitude with high statistical 
accuracy. 

The Vega-referred magnitude system is used throughout 
the paper. 
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2. Observation and Data Analysis 

We performed JHK s -bands imaging observations of 
a part of the Great Observatories Origins Deep Survey 
North (GOODS-N) field with Multi-Object InfraRed 
Camera and Spectrograph (MOIRCS, Ichikawa et al. 
2006) on the Subaru Telescope on 2006 April 4-9,18 and 
May 10-11,18-19. MOIRCS has the field of view of about 
4x7 arcmin 2 with 0.117 arcsec pixel scale. We observed 
one field of view of MOIRCS centered on 12 ?i 36 m 46 s .62, 
+62°13T5".6 (J2000), which included the original Hubble 
Deep Field North (HDF-N, Williams et al. 1996). Here 
we use the high-quality J and iG-bands data sets which 
include only the frames with the seeing size smaller than 
0.5 arcsec (FWHM). The corresponding exposure times 
are 5.0 hours in J-band and 7.7 hours in if s -band. 

The data were reduced in a standard way using the 
IRAF software package. At first the self-skyflat frame was 
made and used for the flat-fielding. We then performed 
the sky subtraction, and co-registered and combined the 
data. The details of the data processing and the data 
quality are described in Konishi et al. (in preparation). 
The FWHMs of the PSFs of the combined images are 
0.42 arcsec in J-band and 0.40 arcsec in if s -band, and 
the if s -band image was convolved with a Gaussian kernel 
to match the PSF to that in J-band. FS23 and FS27 
in the UKIRT faint standard stars were used for the flux 
calibrations in the J and if -bands. We excluded the edges 
of the images with lower sensitivities, and our final multi- 
band combined images have a slightly smaller field of view 
of about 24.3 arcmin 2 . 

Source detection was performed in the X s -band im- 
age using the SExtractor image analysis package (Bertin 
& Arnouts 1996). We adopted MAG_AUTO from the 
SExtractor as the total if-band magnitudes of the de- 
tected objects. For the color measurements, we used the 
fixed aperture with 0.85 arcsec diameter (2xFWHM), and 
the same apertures were used for J and if -bands. 

The 5cr limiting magnitude at 0.85 arcsec diameter is 
about J = 24.6 and K — 23.2, respectively. We estimated 
the la background fluctuation in each band by directly 
measuring the sky fluxes with 0.85 arcsec apertures ran- 
domly placed on the images. The estimated background 
fluctuations were used to calculate the photometric errors 
for the detected objects. We performed the simulations 
using the IRAF/ ARTD ATA package to quantify the de- 
tection completeness in X-band. Our source detection is 
nearly complete for the point source to K ~ 23, and the 
completeness is 90% at K ~ 23.3. We also tested sensitiv- 
ity to the false detections by running SExtractor on the 
inverted if s -band image. Only 15 spurious objects were 
extracted at K < 23, while 1596 objects with K < 23 were 
detected in the normal image. 

3. Results 

3.1. Selection of DRGs 

Figure 1 shows J—KvsK colour-magnitude diagram of 
the objects in the MOIRCS deep imaging field. We have 
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Fig. 1. Color-magnitude diagram of the objects in the 
MOIRCS deep imaging field in the GOODS-North. Galaxies 
with J — K > 2.3 are plotted by large circles. Vertical dot- 
ted-dash line represents the completeness limit at _R"-band. 
Dashed line shows the 2a limit at J-band. For objects with 
J-band fluxes lower than the limit, the lower limits of J — K 
color are plotted by symbols with arrows. Those objects do 
not necessarily lie on the dashed line because the aperture 
sizes for the J — K color and the A"-band magnitude are dif- 
ferent (0.85 arcsec diameter and the Kron aperture, respec- 
tively). Squares indicates the objects which are identified 
as point sources in the publicly available GOODS HST/ACS 
2-band data. 

selected DRGs by the criterion of J — K > 2.3. These 
DRGs are plotted as the large circles, while the others 
are showed as the small crosses. The squares show the 
objects which were identified as point sources in the pub- 
licly available GOODS HST/ACS z-band data (Giavalisco 
et al. 2004). The sequence of these objects is seen at 
J — K < 0.9, which is consistent with the expected locus 
of the galactic stars. This suggests that the flux calibra- 
tion of our data is good enough for the selection of DRGs. 
The 2(7 limit at J-band and the completeness limit at K- 
band are also showed (dashed line and dotted-dash line). 
The deep J-band data (the 2a limit of J = 25.6) allows us 
to sample the objects with J — K > 2.3 down to K ~ 23.3. 

91 (86 at K < 23) DRGs were detected in our field of 
~24.3 arcmin 2 . The cumulative surface density of DRGs 
down to K =21, 22, 23 is tabulated in Table 1. The errors 
of the surface densities are based on Poissonian. The sur- 
face density of 1.15±0.22 arcmin -2 at K < 21 is consistent 
with those in the HDF-South (HDF-S) and MS1054-03 
field reported by Forster Schreiber et al. (2004). 

3.2. Number Counts of DRGs 

Figure 2 shows the number counts of the DRGs in the 
MOIRCS deep imaging field. We show the number counts 
of all if-selected objects and the DRGs in the top panel. 
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Fig. 2. (a) Number counts of all i^-selected objects (squares) and DRGs (circles) in the MOIRCS deep imaging field. The detection 
completeness in iC-band is also shown (dotted line). The vertical dotted-dash line shows our completeness limit. ii"-band number 
counts in other fields are also shown (HDF-South from Labbe et al. 2003, MS1054-03 field from Forster Schreiber et al. 2006, 
GOODS-South from Grazian et al. 2006, Subaru Deep Field from Maihara et al. 2001, Subaru Super Deep Field from Minowa et al. 
2005, high Galactic latitude fields from Bershady et al. 1998). (b) Comparison of DRG number counts. The DRG number counts 
in other fields are also shown (HDF-South, MS1054-03, GOODS-South). The dashed line shows the linear fit for the data points at 
K < 22 in all four fields (this study, HDF-S, MS1054-03, GOODS-S). Errorbars are based on Poissonian (i.e., square roots of the 
observed numbers). Dotted line shows the detection completeness for extended objects with FWHM=0.58 arcsec (see text). 
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Table 1. Number density of DRGs in the MOIRCS deep field 



K magnitude 


Number 


density (arcmin -2 mag 1 ) 


< 18.5 








18.5-19.0 


2 


0.16 ± 0.12 


19.0-19.5 


3 


0.25 ± 0.14 


19.5-20.0 


4 


0.33 ± 0.16 


20.0-20.5 


7 


0.58 ± 0.22 


20.5-21.0 


12 


1.07 ± 0.30 


21.0-21.5 


5 


0.41 ± 0.18 


21.5-22.0 


24 


1.98 ± 0.40 


22.0-22.5 


18 


1.48 ± 0.35 


22.5-23.0 


11 


0.91 ± 0.27 






cumulative density 


A-band limit 


Number 


(arcmin -2 ) 


< 21.0 


28 


1.15 ± 0.22 


< 22.0 


57 


2.35 ± 0.31 


< 23.0 


86 


3.54 ± 0.38 



The A-band detection completeness mentioned in the pre- 
vious section is also plotted. Our data are nearly complete 
down to A = 23, and the number counts in Figure 2 are not 
corrected for the detection completeness. We also show 
A-band number counts in other general fields from the 
literature. The number counts of all A-selected objects in 
our data are consistent with those in other fields down to 
K = 23. 

The number of the detected DRGs and their surface 
density in each A-band magnitude bin are also showed 
in Table 1. The interesting result in Figure 2 is that the 
number counts of DRGs turn over at A ~ 22, while the 
A-band number counts continue to increase to at least 
K = 23, which corresponds to the completeness limit. 

In the bottom panel of Figure 2, we compare the number 
counts of DRGs in our data with those in other fields 
(HDF-S, MS1054-03, GOODS-South). The HDF-S data 
by the FIRES survey (Labbe et al. 2003) have the smaller 
area (~4.5 arcmin 2 ) but reach to the deeper completeness 
limit (K ~23.8) than our data. The areas of the GOODS- 
S and MS 1054-03 fields are comparable with or wider 
than this study (~130 and ~25 arcmin 2 , respectively), 
but these data reach to the shallower limit of A ~21.5-22. 

Though these independent four fields have different ar- 
eas and depths, the number counts of DRGs in these fields 
agree well at K < 22. At K > 22, only the HDF-S data 
have the sufficient depth and can be used for the compar- 
ison with out data. The number counts of DRGs in the 
HDF-S do not show the decrease at A > 22, but are still 
consistent with our results within the uncertainty (see also 
the next section). 

In order to demonstrate the turnover of the DRG counts 
or the deficit of the faint DRGs relative to the brighter 
ones, we performed the linear fitting for the number counts 
of DRGs at 18.5 < K < 22 using the data in all four fields. 
The linear fitting at 18.5 < K < 22 seems to be reasonable 
because Foucaud et al. (2006) reported the break feature 
of the DRG counts lies at A ~ 18 and the data over 18.5 < 
K < 22 of all four fields were fitted well with the linear 



line. The result is plotted as the dashed line in Figure 2. 
We also tried to fit the MOIRCS data at K < 22 with the 
Maximum Likelihood method without binning the data, 
and confirmed the very similar result of the fitted slope. 
Thus the uncertainty due to the binning of the data does 
not affect the fitting result, even if there is a bin with the 
relatively low counts (A =21-21.5 bin). 

When we extrapolate this linear line to K ~ 23, the 
surface density of the faint DRGs with 22 < A < 23 in our 
field is clearly deficient. The density in the 22.5 < A < 23 
bin is about a factor of three lower than the extrapolation 
from the number counts at the brighter magnitude. If we 
perform the similar linear fitting for the MOIRCS data 
over 18.5 < K < 23 with the Maximum Likelihood method, 
the fitted slope become 0.19 ± 0.04, which is significantly 
flatter than the value for the fitting at 18.5 < K < 22 (0.33 
± 0.07). It is seen that even in the HDF-S, the number of 
the faint DRGs is lower than the extrapolation from the 
number counts at K < 22. 

4. Discussion 

In this section, we discuss the turnover of the number 
counts of DRGs at K ~ 22 found in the previous section. 

At first, we consider the possible spurious effects which 
could cause the deficit of the faint DRGs. One possibility 
is that the MAG_AUTO would systematically underes- 
timate the total flux in the A-band near the detection 
limit (e.g., Labbe et al. 2003), and this could cause the 
underestimation of the number of the faint DRGs. The 
number counts of DRGs, however, begin to decrease at 
K ~ 22, which is one magnitude brighter than the com- 
pleteness limit in A-band, and this does not seem to be 
the case. For example, Labbe et al. (2003) showed such a 
underestimation occurs at K ~ 24, which is about 0.3 mag 
fainter than the completeness limit in their FIRES data. 
Our simulation also indicates that the significant under- 
estimation of the A-band flux does not exist down to at 
least A = 23 (Konishi et al. in preparation). In fact, the 
A-band number counts continues to increase to A = 23 in 
Figure 2. We also confirmed that the size distribution of 
DRGs with K ~ 23 is not significantly different from that 
of the bluer objects. Therefore there is no reason to think 
that only DRGs are less complete at A ~ 23. 

Since most DRGs in our data are resolved, the com- 
pleteness problems could be more severe than those for 
the point source. Therefore we also checked the detec- 
tion completeness for extended objects. We fitted the A- 
band surface brightness of relatively bright DRGs with 
K < 21 with the Moffat function, and found that the 
average FWHM is about 0.58 arcsec. Using this sur- 
face brightness profile (assuming the same profile for the 
fainter DRGs), we performed the similar simulation with 
that for the point source to quantify the detection com- 
pleteness as a function of A-band magnitude. The result 
is showed in the bottom panel of Figure 2 (dotted line). 
The completeness is more than 90% at A < 22.5 and is still 
80% at A = 23. The detection incompleteness in A-band 
cannot explain the observed deficit of the faint DRGs. 
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Fig. 3. Spatial distribution of the DRGs. The size of the symbols is scaled according to apparent magnitudes in K-band (the bigger, 
the brighter). Dashed line shows the region of the original HDF-N. Dotted-dash lines show the field of view of the MOIRCS data. 
MOIRCS achieves the field of view of 4 X 7 arcmin 2 with two HAWAII-2 arrays. The present version of the combined data were 
separately reduced for each array, and the deep portion of the field of view is divided into two areas. Two corners (upper left and 
upper right) are vignetted at the Cassegrain focus of the telescope. 



Second possibility is that the larger error in J — K color 
at the faint magnitude affects the sampling of DRGs. 
Although the depth of our J-band data (the 2a limit of 
J = 25.6) is sufficient to sample the objects with J — K > 
2.3 to K = 23, some fraction of the faint DRGs could be 
missed due to the photometric error. However, we con- 
clude that the error in J — K does not cause the deficit 
of the faint DRGs, since the similar effect of the contam- 
ination from the bluer objects to the DRG sample also 
exists. As discussed in Foucaud et al. (2006), since the 
overall J — K distribution of objects has the peak at the 
color substantially bluer than J — K ~ 2.3 and the objects 
with J — K> 2.3 are relatively rare (Figure 1), it can 
be expected that the number density of DRGs is rather 
boosted at fainter magnitudes. Therefore these photo- 
metric errors do not seem to cause the observed deficit of 
DRGs at K > 22. 

Although the number counts of DRGs in the HDF-S 
do not show the turnover at K ~ 22 as seen in the previ- 
ous section, the difference of the number density of DRGs 
at K > 22 between in the HDF-S and in our MOIRCS 
field can be explained by the ficld-to-ficld variance, when 
the strong clustering of DRGs is taken into account (e.g., 
Grazian ct al. 2006, Quadri et al. 2006). Grazian ct al. 
(2006) pointed out the large discrepancy of the number 
density of DRGs with K < 22 between the HDF-S and 
HDF-N (e.g., Dickinson et al. 2000, Fontana et al. 2000). 
They reported that only two DRGs (one at K < 21, one 
at 21 < K < 22) exist in the HDF-N, while many DRGs 
are found in the HDF-S. This can be seen in our data. In 
Figure 3, we plot the spatial distribution of DRGs in the 
MOIRCS deep field. The region of the original HDF-N is 



also showed, and there are only two DRGs with K < 22 
and one DRG with 22 < K < 23 in this region. On the 
other hand, there are many DRGs over wide range of K- 
band magnitude outside the HDF-N region, and not only 
the bright DRGs but also the faint ones with 22 < K < 23 
appear to show a relatively strong clustering. Such a spa- 
tial distribution of DRGs indicates that the small survey 
areas of the HDFs could introduce large uncertainty and 
the wide area survey is essential for the estimation of the 
number density of these galaxies. We will present the de- 
tailed analysis of the clustering property of these galaxies 
in Ichikawa et al. (in preparation). 

What does the deficit of the faint DRGs at K > 22 
mean ? The previous studies found that most DRGs exist 
at 2 < z < 4 (van Dokkum et al. 2003, Reddy et al. 2005, 
Forster Schrcibcr et al. 2004). Our photometric redshift 
estimation with the HST/ACS and MOIRCS data also 
suggests the similar redshift distribution; 69 out of 86 
DRGs with photometric redshift (those detected in more 
than two bands) have 2 < z p hot < 4 (Konishi et al. in 
preparation). Keeping in mind that some significant frac- 
tion of bright DRGs confirmed to be at z < 2 (e.g., Daddi 
et al. 2004, Conselice et al. 2006) and that the reliability 
of photometric redshift for the faint DRGs has not yet 
confirmed by spectroscopy, we assume that most DRGs in 
our data lie at 2 < z < 4 in the following discussion. 

The observed color of J — K = 2.3 roughly corresponds 

to the rest-frame U - V ~ 0.4 (U - B 0.1) at 2 < z < 4, 

which is similar with those of Sc-Sd galaxies seen in the 
present universe. Such a red color for high-z objects 
(about 2-3 Gyr age of the universe at 2 < z < 4) is con- 
sidered to be produced by Balmer/4000A-break of rela- 
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Fig. 4. A'-band number counts of galaxies with 2 < z p hot < 4 
in the MOIRCS deep imaging field. Circles show those galax- 
ies with J — K > 2.3 (i.e., DRGs) and squares represent those 
with J — K < 2.3. Errorbars are based on Poissonian. The 
completeness limit is also shown (dottcd-dash line). 

tively old stellar population and/or heavy dust extinction 
as mentioned in Section 1. Therefore, the turnover of the 
number counts of DRGs indicates that while there are 
many bright (massive) galaxies with the old population 
and/or dusty star formation activity, the number of the 
faint galaxies with the similar red color is relatively small 
at that epoch. Labbe et al. (2005) investigated the stel- 
lar mass and SED of DRGs and LBGs with K < 22 and 
found that DRGs dominate the high-mass end at z ~ 3. 
On the other hand, our result from the wide and deep 
data could suggest that the low-mass end is dominated by 
galaxies with blue colors, which do not satisfy the criterion 
of DRG. In Figure 4, we show the if-band number counts 
of galaxies with 2 < z p hot < 4 for DRGs ( J — K > 2.3) and 
those with J — K < 2.3, separately. It is seen that the frac- 
tion of those with J — K < 2.3 increases with decreasing 
.ftT-band flux at K > 22, while the number density of DRGs 
is comparable to that of the bluer objects at K ~20-21. 
Kajisawa and Yamada (2005, 2006) showed that most low- 
mass galaxies have the blue rest- frame U — V color and the 
more massive galaxies tend to have the redder color even 
at z ~ 2.5. Kajisawa & Yamada (2006) also evaluated 
that the transition mass between these red and blue pop- 
ulations is about 6x10 10 Mq at z ~ 2.5, and in fact most 
of galaxies with M* < 1O 1O M0 have the rather blue color 
of U — V ^5 0-1 (their Figure 2). Since the magnitude of 
K = 22 of galaxies with J — K ~ 2.3 corresponds to about 
l-6xlO lo M0 at 2 < z < 4 (we calculated the stellar mass 
as in Kajisawa & Yamada 2006, using the GALAXEV 
population synthesis model of Bruzual & Chariot 2003), 
the deficit of DRGs at K > 22 could be explained by such 



a mass-dependent color distribution, if the similar trend 
continues to z ~ 4. 

We plan to perform the statistical analysis of the SEDs 
of these galaxies over the wide range of luminosity (and 
mass) in order to constrain their physical properties in a 
forthcoming paper. 
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